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By combining neutron inelastic scattering (NIS) and first-principles calculations, we have investi- 
gated the lattice dynamics of MOFS. The structural stability of MOFS was evaluated by calculating 
the three cubic elastic constants. We find that the shear modulus, C44 = 1.16 GPA, is unusually 
small, while two other moduli are relatively large (i.e. cn = 29.42 GPa and C12 = 12. S6 GPa). We 
predict that MOFS is very close to structural instability and may yield interesting new phases under 
high pressure and strain. The phonon dispersion curves and phonon density of states were directly 
calculated and our simulated NIS spectrum agrees very well with our experimental data. Several 
interesting phonon modes are discussed, including the softest twisting modes of the organic linker. 

PACS numbers; 63.20.-e, 62.20.Dc, 61.82.Pv 



X3 

o 



> 

(N 

o\ 

O 

o 



I 

O 

o 



Metal-organic framework (MOF) compounds fiiSi^i^ 
which consist of metal-oxide clusters connected by or- 
ganic linkers (see Fig. 1), are a new class of nanoporous 
materials with very promising potential applications such 
as energy storage, gas separation, and template synthesis 
of nanoclusters and catalysts. There are an almost expo- 
nentially growing number of studies on MOF materials, 
mostly focusing on the optimization of the metal-oxide 
clusters and/or the organic linkers to improve the gas 
adsorption properties. Much less attention is paid to the 
fundamental properties of this interesting class of materi- 
als such as their structural stability and lattice dynamics. 
Understanding the stability and dynamical properties is 
clearly needed in order to optimize these materials for de- 
sired properties. For example, it was recently proposed 
that the MOF structure can encapsulate Ceo molecules 
and result in superconductivity upon doping"". In such 
hybrid structures, it is important to know the phonon 
spectrum of the MOF lattice and its coupling to the 
electronic structure. Furthermore, a lack of knowledge 
concerning MOFs lattice dynamics presents a significant 
obstacle in the quantum dynamics study of gas molecules 
(H2, CH4 etc.) adsorbed on MOFs&i 

Here we report a detailed study of the structural stabil- 
ity and lattice dynamics of MOFS (the most widely stud- 
ied MOF material) from combined neutron inelastic scat- 
tering (NIS) and first-principles calculations. The struc- 
ture of M0F5 (Fig.l) is highly symmetric (space group 
Fm-3m) and consists of Zn04 clusters at the tetrahedral 
site of an fee lattice, linked by 1,4-benzenedicarboxylate 
(BDC) to form the three dimensional framework. 

Our first-principles calculations were performed within 
the plane-wave implementation of the local density ap- 
proximation to density functional theory in the PWscf 
package^. We used Vanderbilt-type ultrasoft potentials 
with Perdew-Zunger exchange correlation. A cutoff en- 
ergy of 408 eV and a 1x1x1 k-point mesh were found 
to be enough for the total energy to converge within 
0.5meV/atom. 




FIG. 1: (Color online) The crystal stucture of MOFS, which 
consists of BDC-linkers connecting Zn04 clusters located at 
the tetrahedral site of an fee lattice. For clarity, some of the 
atoms were not shown. 



We first describe the dependence of the structure on 
the cell volume (i.e. external pressure). Figure 2 shows 
how the lattice parameter and the total energy of the 
cell changes with external pressure and the total volume, 
respectively. We optimized the atomic positions for each 
volume and did not find any structural instability over 
the pressure range studied. From -0.5 to 0.5 GPa, we find 
a linear pressure dependence of the lattice constant with 
a slope of 0.46 A/GPa. Fitting the energy vs. volume 
curve to the Murnaghan equation of stated yields a bulk 
modulus of 18.2 GPa. The minimum energy corresponds 
to lattice constant ac — 25.58 A, which is in excellent 
agreement with the experimental valueiS of a = 25.91 A. 
The optimized atomic positions were also found to agree 
well with the experimental values. 

A better insight into the structural stability can be 
gained from the elastic constants, which we calculated 
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FIG. 2: (Color online) Left: The cubic lattice parameter 
of M0F5 versus pressure, indicating linear dependence near 
equilibrium structure. Right: The total energy versus primi- 
tive cell volume (dots) and the fit by Murnaghan equation of 
state (solid line). 
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FIG. 3: (Color online) The three cubic elastic constants of 
MOF5. The dots are the actual calculations and the solid 
lines are the quadratic fit. The deformation matrices for each 
distortion are also shown. The right-bottom panel summa- 
rizes the calculated elastic moduli and the sound-velocities. 
Here p is the mass density of MOF5 (0.59 g/cm^). 



following the standard scheme for a cubic crystalii*i^, 
starting from the fully optimized M0F5 structure. Three 
types of strains were applied to the conventional fee cell of 
M0F5 (see Fig. 3). Strain / is a stretch along the z axis; 
strain // is a monoclinic shear about the z axis. Both 
strain / and // are volume-conserving. Strain /// simply 
changes the lattice constant a. For all three strains, there 
is only one variable, 7, in the strain tensor, reflecting the 
magnitude of the strain. All calculations were done on 
the strained primitive cell, which is structurally equiva- 
lent to the strained conventional cell. Atomic positions 
were fully optimized and E vs. 7 data were obtained. 



E(^)/V was then fit to a polynomial in 7, where V is the 
unit cell volume. The results are summarized in Fig. 3. 
As expected, the calculated elastic constants satisfy the 
stability criteria (i.e. cn > C12, cn -|- 2ci2 > 0, and C44 > 
for a cubic crystal). However, the shear modulus C44 
is very small (1.16 GPa), indicating that M0F5 is very 
close to structural instability. Besides, cn — C44 is sig- 
nificantly larger than zero and according to the Cauchy 
relationi^, this indicates a significant deviation from a 
central intermolecular potential. The small C44 also sug- 
gests that M0F5 could collapse into a potentially useful 
structure under certain shear-stress. For examples, the 
pristine M0F5 structure has a large pore volume but the 
pores are too open to hold hydrogen molecules at high 
temperatures; a collapsed M0F5 structure could have 
finer cavities connected by smaller channels, which are in 
principle better for H2 storage. 

After having shown that the M0F5 structure is sta- 
ble with a very small shear modulus, we now discuss 
the lattice dynamics properties at the zone center and 
along high-symmetry directions in the Brillouin zone. 
The phonon density of states (DOS) and dispersion 
curves were calculated using the supercell method with 
finite differenceii^ The primitive cell was used and the 
full dynamical matrix was obtained from a total of 26 
symmetry-independent atomic displacements (0.03 A). 

The primitive cell of M0F5 contains two formula units 
(Zn40i3-(C8H4)3) giving rise to a total of 318 phonon 
branches. The phonon modes at F are classified as 

r (q=0) = 9Aig(R) + 3Ai„ + 3^29 + 9^2u + l2Eu 
+ 12Sg(R) + imu + 24T2g(R) + 24Ti„(IR) + ITTig 

where R and IR correspond to Raman- and infrared- 
active, respectively. The crystal symmetry implies 105 
Raman- and 72 IR-active modes. In Table I, we list 
the calculated energies at F. We hope that our cal- 
culations will initiate more experimental work such as 
Raman/IR measurements to confirm the gamma phonon 
energies that we calculated here. In the absent of such 
high-resolution data, we have performed neutron inelas- 
tic scattering measurements of the phonon DOS and com- 
pared it with our calculations. 

The experimental NIS spectrum shown in Fig. 4 was 
collected on a powder sample of M0F5 using the Fil- 
ter Analyzer Neutron Spectrometer-^ (FANS) under con- 
ditions that provide fuU-width-at-half-maximum energy 
resolutions of 2-4.5% of the incident energy over the range 
probed. To compare the NIS data with theory, the NIS 
spectrum was computed for a 10x10x10 q-point grid 
within the incoherent approximation^'* with instru- 
mental resolution taken into account. As shown in Fig. 
4, the agreement between calculation and experiment is 
excellent. In contrast to Raman and IR spectroscopy, 
neutron spectroscopy is not subject to any selection rules, 
thus all modes falling in the energy range covered by our 
measurement (~20 to 170 meV) were detected. Within 
this E range, the NIS spectrum is dominated by one- 
phonon processes. The two-phonon contribution coming 
from the combination of one-phonon processes is not sig- 
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TABLE I: List of phonon symmetries and calculated energies 
(in meV) of M0F5 at the F point of the primitive cell. 
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FIG. 4: (Color online) Measured NIS spectrum of M0F5 
at 4 K (top curve) along with calculated spectrum (bottom 
curves). The calculated 1 (dotted line) and 1+2 phonon con- 
tributions are shown. 

nificant. Also note that the NIS spectrum is dominated 
by hydrogen displacements because of the relatively large 
H neutron scattering cross section. 

Figure 5 shows the dispersion curves of the low-energy 
modes along high-symmetry directions in the Brillouin 
zone. The high energy modes are quite dispersionless 
and thus not shown. It is interesting to point out that 
some of the modes have significant dispersion, such as 
those between 25-30 meV. This is surprising, considering 
the rigidity of the units (i.e., the Zn04 cluster and the 
BDC linker). We find that the dispersion comes from the 
points where these two rigid units are connected. This 
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FIG. 5: Calculated phonon dispersion curves along high- 
symmetry directions in the Brillouin zone for MOF5. For 
clarity, only low energy (0-40 meV) portion is shown. 



is also the reason for the very low shear modulus C44. 
Another interesting result is the very low energy modes 
near 3 meV, which are quite flat. Some of these phonon 
modes will be discussed in detail below. Up to now, no 
phonon-dispersion measurement on M0F5 has been re- 
ported, partly due to the difficulty of growing a large 
single crystal needed for the experiment. We hope that 
in the near future, some of the phonon dispersions that 
we predict here can be experimentally confirmed. 

We next discuss several interesting modes revealed by 
the phonon calculations. Animations of these phonon 
modes are available in Refiii for direct visualization. The 
lowest energy mode {A2g, 2.59 meV) is singly degenerate, 
with all C6II4 linkers twisting around the crystal axis di- 
rection. This mode is basically dispersionless (see Fig. 
5) which is due to negligible interaction between two ad- 
jacent benzene linkers along the crystal axis. The second 
lowest energy mode (Tig, 3.12 meV, 3-fold degenerate) 
can be described as two of three C6H4 linkers twisting 
around the two corresponding crystal axes. These two 
softest twisting modes imply that the bonds between the 
organic linkers and the metal clusters are the most unsta- 
ble part of the MOF structure. The ^^0.53 meV energy 
difference between the two modes suggests a coupling 
of the twisting motions of the nearest-neighbor benzene 
linkers. 

Since the energy of the twisting mode is very low, one 
may wonder if the benzene ring is free to rotate at high 
temperatures. This could have important consequences 
on the hydrogen absorption dynamics in M0F5. If the 
linker is rotated 90°, it could block the channels connect- 
ing the cubic cavities in M0F5 available for hydrogen 
storage. Motivated by these arguments, we have investi- 
gated the dynamics of the linker using quasielastic neu- 
tron scattering measurements^^. However, we did not 
observe any free rotation of the linker at high tempera- 
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FIG. 6: (Color online) The variation of the total energy as 
one of the C6H4 groups (part of the BDC organic linker of 
MOF5) is twisted around the crystal axis direction. Note 
that the initial twisting does not cost much energy (see inset) 
although the energy barrier for the 180° flipping (two-site 
jumping) of the C6H4 group is rather large (~537 meV). 

ture nor any evidence of the 7r-rotational jumps of the 
linker. These findings later became clear when we looked 
at the potential curve as the C6H4 linker is rotated. We 
performed total energy calculations as a function of the 
rotational angle of the benzene plane. For simplicity, we 
only twisted one of the six benzene rings in the prim- 
itive cell. The results are shown in Fig. 6. The over- 
all barrier is surprisingly large (~537 meV), consistent 
with the quasielastic neutron scattering results. How- 
ever, twisting for a few degrees only costs an energy of 
several meV (see the inset of Fig. 6), in agreement with 
the phonon modes being very soft. The third lowest en- 
ergy mode (3.50 meV, 3-fold degenerate) involves 
all C6H4 linkers translating along crystal axis a, 6, and 
c. All other modes with E < 23 meV are various com- 
binations of the motion of the C6H4 linker as a whole 
and the vibration of the Zn04 clusters. Two modes at 
slightly higher energy (23.69 meV, and 28.67 meV, 



Aig) are the breathing modes of the Zn04 cluster, with 
the clusters vibrating in the unit cell out-of-phase and 
in-phase, respectively. Most of other modes falling in 
the energy range of our NIS measurement are one of the 
twisting, stretching, breathing, bending and wagging mo- 
tions of the benzene rings. They can be roughly divided 
into two groups. With E < ~70 meV, the phonon modes 
are more dominated by the C displacements of the C6H4 
linkers. The modes in the higher energy region are more 
dominated by H displacements. 

The highest energy modes (380-385 meV, not covered 
by our NIS data) are totally dominated by the C-H cova- 
lent bond and all are associated with various C-H stretch- 
ing modes of the C6H4 group. 

In summary, we presented a detailed lattice dynam- 
ics study of M0F5 combining first-principles calculations 
and neutron inelastic scattering measurements. Our re- 
sults indicated that MOFs are an interesting class of 
materials in terms of structural stability and dynami- 
cal properties. We identified the modes associated with 
the twisting of the linker to be the softest modes in the 
structure and therefore expected to be populated first 
with increasing temperature, eventually leading to struc- 
tural instability. This is also consistent with the very 
small shear modulus C44 that we have calculated. De- 
spite the very low energy of the linker twisting modes, 
we found a very large energy barrier for the full rotation 
of the linker. We hope that our calculations will initiate 
further experiments to confirm our predictions. In par- 
ticular, our results suggest that M0F5 is very close to 
structural instability, possibly leading to an interesting 
structural transformation under high pressure. 
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